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ABSTRACT 

We investigate the muon neutrino event rate in km^ neutrino telescopes 
due to a number of galactic supernova remnants expected on the basis of these 
objects' known 7-ray signals. We evaluate the potential of these neutrino signals 
to exhibit evidence of the sub-dominant neutrino oscillations expected in various 
neutrino mixing schemes including pseudo-Dirac scenarios and the Exact Parity 
Model. With ten years' data, neutrino signals from Sgr A East should either 
discover or exclude neutrino oscillations governed by a Sm^ parameter in the 
range 10~^^ to 10~^^ eV^. Such a capability is not available to terrestrial or solar 
system neutrino experiments. 

Subject headings: acceleration of particles — cosmic-rays — elementary particles: 
neutrinos — radiation mechanisms: nonthermal — supernova remnants 
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1. Motivation and Plan 

In this work we present a novel extension of the basic idea that astrophysical neutrinos 
can probe tiny values of the difference in the squared masses of relevant neutrino mass 
eigenstates, Sm^. This, namely, is that through the observation of a deviation away from 
pure power-law scaling - in other words, a spectral distortion - in a particular, galactic 
supernova remnant's observed muon neutrino spectrum, an experimentalist can infer the 
existence of exactly such a tiny mass splitting. Such mass splittings are generic to a number 
of extensions of the Standard Model such as the various active <-> sterile pseudo-Dirac 
scenarios and the Exact Parity Model. The crucial advance in this proposed method is that 
one does not need to observe the species of neutrino into which the neutrinos emitted at 
the source are oscillating in order to diagnose oscillations. This is important because both 
astrophysical i/gS and ly^s are expected to be considerably more difficult to detect than 
astrophysical i^^u's, at least at more moderate energies. 

The plan of the paper is as follows: in §2 we discuss recent developments, both 
empirical and theoretical, concerning galactic supernova remnants (SNRs), including recent 
considerations of the factors limiting the maximum energies to which they may accelerate 
particles. In §3 we discuss neutrino production at SNRs through pion decay and the 
relationship between an SNR's 7-ray signal and its expected neutrino flux at Earth. In §4 
we briefly review neutrino oscillations and the status of the various experiments purporting 
to demonstrate such oscillations. We also introduce here the idea of tiny 5m^'s and review 
their theoretical motivations. §5 considers the effect such tiny 5m^'s might have on the 
phenomenology of SNR neutrinos. In §6 we briefly review the extensive code we have 
written to model neutrino telescope detection of SNR neutrinos. Finally, we set out the 
results of this code in §7 which goes on to demonstrates the feasibility of employing the 
spectral distortion method briefly described above (and set out in further detail later) to 
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search for very long wavelength neutrino oscillations. 



2. Particle Acceleration in SNR Shells 

Cosmic-ray ions and electrons up to (and possibly exceeding) energies of ~ 10^^ eV 
(near the so-called 'knee' in the distribution observed at earth) are widely believed to 
be produced by galactic SNRs. The dominant acceleration mechanism in SNRs appears 
to be diffusive (or first-order Fermi) acceleration at the remnants' forward shocks. Some 
support for this supposition was provided by the EGRET experiment aboard the Compton 
Gamma-Ray Observatory^ which detected a large number of (at first) unidentified sources 
in the super-50 MeV band, both in and above the Galactic plane. Six of these EGRET 
sources turned out to have compelling associations with relatively young SNRs QEsposito 



et al. 1996|) , with gamma-ray fluxes at earth of typically a few x 10 photons cm ^ s ^ 



above 100 MeV. 

The inference to be drawn here is that the accelerated ion (mostly proton) and 
electron distributions of the EGRET SNRs spawn neutral and charged pion decay and the 
accompanying bremsstrahlung, inverse Compton, and synchrotron emission, that account 
for their broad-band photon spectra from radio to gamma-ray energies. In particular, the 
high energy photon spectra of these objects inferred from the EGRET measurements show 
deviations from pure power-law scaling around 100 MeV that have been argued to be 
evidence of vr" decay ( [Gaisser, Protheroe, fc Stanev 1998| ; [Markofl', Melia fc Sarcevic 1997| ). 



These vr^'s would be expected from collisions between shocked p's in the SNRs' expanding 
shells and ambient nucleons. 

EGRET has also detected similar 7-rays from the Galactic Center (GC; |Mayer-| 
Hasselwander et al. 1998) that Melia and others have strongly argued are associated with 
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TT 



decay at the SNR-like object Sgr A East ([Melia et al. 1998 ; Markotf, Melia fc Sarcevi" 



1997|) . (See Crocker, Melia and Volkas 2000 for an analysis of the significance of Sgr A East 



for V astronomy.) 

The problem with this otherwise straightforward interpretation, however, is that a 
simple extension of the EGRET SNRs' GeV spectra to higher energies produces a TeV flux 
in excess of the upper limits estabhshed by the Whipple ( Puckley et al. 1998| ; [Rowell et] 
al. 2000|) and High-Energy Gamma-Ray Array (HEGRA; [Prosch et al. 1995|) atmospheric 



Cerenkov telescopes and the CYGNUS extensive air shower array ([Allen et al. 1995| ; except 
in the case of the GC which has been seen - marginally - by Whipple; [Buckley et al. 19IF7| ). 
Thus, a crucial question arises as to whether shell-type remnants can indeed supply the 
observed Galactic cosmic-ray population, while at the same time explaining the emission 
from the handful of unidentified, EGRET sources associated with SNRs. It may well be 
that the physical conditions conducive to producing large numbers of relativistic particles 
also work against energizing any one of these to TeV, let alone PeV, energies. 

In their detailed analysis of this phenomenon. Baring et al. (1999) adopted the view 
that the maximum energy attainable by ions in diffusive shock acceleration is determined by 
two complementary conditions, viz: (1) that the acceleration time (as a function of energy) 
ought not to exceed the age of the remnant (for the free expansion or early Sedov phase), 
and (2) that the diffusion length of the highest energy particles ought not to exceed some 
fraction of the shock radius. They concluded that, for SNRs in a homogeneous environment, 
the circumstances that favor intense gamma-ray production in the EGRET and sub- TeV 
bands (namely an ISM density > 1 cm~^) limit acceleration of particles to energies well 
below the cosmic-ray knee. Thus, dense remnant environments, as might be expected for 
the low Galactic-latitude sources in the Whipple and HEGRA surveys, produce luminous 
emission in the GeV range with a simultaneous absence of TeV 7-rays because the resultant 
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spectrum cuts-off (or significantly steepens) before tliis energy. 

At tlie same time, tlie detection of TeV emission from sources out of the Galactic 
plane (and hence presumably in regions of lower ISM density) leaves little doubt that 
there exists somewhat of an anti-correlation between GeV gamma-ray luminosity and 
super- 10 TeV cosmic-ray production in individual sources. Neither the remnant of SN 1006 
(Tanimori et al. 1998), nor SNR RX J1713.7-3946 (Muraishi et al. 2000), were detected 
by EGRET, yet their gamma ray signals were observed by CANGAROO. (We present the 
high energy 7-ray data for 7-Cygni and SN 1006 in Figure (1) for comparison.) Analysis of 
the TeV gamma-ray emission from the northeast rim of SN 1006 indicates either super-50 
TeV electrons inverse Compton scattering 2.7 K cosmic microwave background photons 
or the decay of 7r° mesons - particularly if the remnant is close by ([Aharonian fc Atoyan] 
1999|) - or both. Note that the putative, high-energy electrons might be either directly 
shock accelerated or decay products resulting from collisions of shock-accelerated protons 
and nuclei. Given that the 7-ray emission is from the rim, they may not, however, be 
accelerated by electro-magnetic processes associated with the remnant's neutron star. In 
any case - whether the gamma-rays are hadronic or leptonic in origin - SN 1006 (and 
SNR RXJ1713. 7-3946 ) present evidence for shock acceleration of protons to high energies 
(at least ~ 100 TeV) given that more significant loss processes act on shock-accelerated 
electrons than protons in SNR environments. 

It appears, therefore, that should the cosmic-ray distribution be produced by SNR 
shells, the most likely sources for the ~ PeV particles, like SN 1006, lie in regions 
of relatively low ISM density. For SN 1006 the TeV luminosity determined from the 
CANGAROO observations is of order lO^'^ erg.s"^, the exact figure depending on distance. 
Such energies are highly desirable for the neutrino physics we wish to probe in this analysis. 
We go on to calculate the level of neutrino emission from relevant SNR sources, which we 
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expect will be predominantly out of the Galactic plane. 

3. Neutrino Production at Supernova Remnants 

Any proton-nucleon collision process at energies sufficient to generate pions leads to 
both z/g and z/^ production from vr^ and /i^ decay; tt^ — > /i^ and /x^ Ve i^fj.- Note 

that we take u to mean u ot V here, as we shall often do in the remainder of this paper, 
because neutrino telescopes will not be able to distinguish between u and V at the energies 
attained by SNR neutrinos. The ratio of neutrino flavors near the point of generation in an 
SNR will be very close to z/g : z/^ : z^^ = 1 : 2 : 0, as simple channel counting in the above 
reactions would indicate. Certainly we expect no significant Ur component at the SNRs 
(from charmed hadron decay) given the energies attained by the shocked protons (almost 
certainly no more than 10 PeV). 

Given that these pion production and decay processes are well-understood physics, 
the high-energy 7-ray signals from an SNR can, in fact, provide an empirical handle on 
its expected total neutrino emission, Q,j{Ey). In more detail, Qij{Ey) can be expressed as 
a function of the 7-ray emission at the SNR, Q^{E^), the index of the SNR's power-law 
proton spectrum above some energy (a ~ 2 such as would result from shock acceleration), 
and r = {m^/mT^y ( |Blasi and Melia 2000| ). Note here that the decay neutrino spectrum 
will have the same spectral index, a, well above the A resonance domain, i.e., for energies 
^2 — 3 GeV. We effect the relative normalization between the 7-ray and neutrino 
distributions at 10 GeV where the pions take on a power-law form. From this we can derive 
a z/^ flux at earth and consequent event rate in a km^ telescope (piasi and Melia 2000 ). 

As a preliminary exercise then, one may determine whether the muon neutrino signal 
from a particular SNR is above atmospheric neutrino background at an energy less than 
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the maximum energy of the SNR i/'s. This is about 1/12 the maximum energy to which the 
SNR might shock p's (~ 10 PeV) given that {xp)-, the average momentum carried by the 
secondary pions relative to the parent proton is ^ 0.2 ( [Alvarez- Mufiiz et al. 2000| ; though 
note that we shall find, in fact, given the expected power law distributions of neutrino 
energies, that low statistics prevent observation of neutrino signals up to the maximum 
cut-off anyway). To perform this preliminary calculation, we assume a neutrino spectrum of 
spectral index a = 2 - a value well-justified theoretically for the case of shock acceleration 
in SNR environments^ and, in many cases able to be inferred directly from an SNR's photon 
spectrum - normalized to its measured or inferred 10 GeV 7-ray signal. We also assume 
a reasonable value of ~ 1° for a neutrino telescope's angular resolution. Granted these 
inputs, we find that a window on the neutrino signal from these objects is opened and we 
can consider what particle and astro- physics we might extract from this ultra-high energy 
and ultra-long baseline neutrino signal. 

The calculations we have performed show that SNRs typically emerge above the 
atmospheric neutrino background at around a TeV, normalizing to their 10 GeV 7-ray 
fiuxes and assuming a 1 km'^ detector. Given that this is at or above the energy regime at 
which the six EGRET SNRs (but not Sgr A East) can be expected to cut-out, we cannot 
safely predict significant neutrino signal from these objects in a km^ detector. On the other 
hand, the two CANGAROO sources can be expected to generate around 5 muon-like events 
per year in such a detector and Sgr A East 50 muon-like events, where we have factored in a 
fiux attenuation of one half due to the averaged z/^ Vr oscillations we expect on the basis 
of the SuperKamiokande atmospheric neutrino data ( [Fukuda et al. 1998a,b,c|) - see later. 



^Furthermore, such a value is concordant with that expected for the parent spectrum of 
galactic cosmic-rays as determined by standard 'leaky box' models ( [Aharonian fc Atoyan 
19991) . 
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We now go on to consider the issue of neutrino oscillations in somewhat greater detail. 



4. Neutrino Oscillations Between SNR and Earth 

4.1. Distance Considerations 

The distance between the neutrino source and detector is about 8 kpc ~ 2.5 x 10^^ cm 
in the case of the Sgr A East. Distance determinations to SN 1006 lie between around 0.7 
kpc ( IWillingale fc West 1996| ) to 2 kpc ( [Winkler fc Long i99l\ ). For SNR RXJ1713. 7-3946 



they likewise range from 6 kpc ( |Slane et al. 1999 ) to 1 kpc ( Koyama et al. 1997 ) 



We take the GC neutrino source - Sgr A East - to have a linear dimension of around 
10 pc ~ 3 X 10^^ m. Given the range of distance determinations, the linear dimensions of 
SN 1006 and SNR RXJ1713. 7-3946 range between 6 and 27 pc and 7 and 41 pc respectively. 
These distances are relevant because we need to know how neutrino oscillation lengths 
compare with the size of the emitting object to determine whether the neutrino source is 
flavor coherent. If the former are small compared to the latter, then, because neutrinos 
are emitted from points distributed across the source, the oscillations will be averaged out. 
Alternatively, if the latter are large compared to the former, then no averaging due to 
the flnite size of the source will be needed and the source is essentially flavor coherent for 
neutrinos of a given energy. 

Note that two types of averaging generally need to be done: over distance, and over 
energy. We need only consider distance averaging due to the flnite size of the u source. 
For SNRs the source distance scales involved are at least six orders of magnitude larger 
than those for the detector (1 A.U. ~ 5 x 10~^ pc). i^etector-dependent distance averaging, 
then, will not impact on calculations concerning neutrinos from SNRs. We do not address 
the issue of energy averaging due to the flnite energy resolution of the detector in great 
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detail in this paper (though we shall only ever consider the binning of event data into the 
coarse-grained scale of a decade of energy - see later) . 



4.2. Introduction to Neutrino Oscillations 

For the moment we consider 2-flavor oscillation modes ^ for illustrative purposes. 
Suppose a beam of flavor a is produced at x = 0. Then at a point x distant from the source 
the oscillation probability is 



P(a -^13) = sin^20sin2 f^?— ) 



The parameter 9 is the 'mixing angle' which determines the amplitude of the oscillations. 
The value 6 = 7r/4, which leads to the largest possible amplitude, is termed 'maximal 
mixing'. The parameter Lose is the 'oscillation length' - the length scale over which 
oscillation between two mass eigenstates occur - and is given by Lose = This works 

out to be 

E/PeV 

Note that the oscillation length increases linearly with energy. The parameter 

(5m^ = \ml — m^l is the squared-mass difference between the two mass eigenstate neutrinos. 

Totally averaged oscillations see the x-dependent sin^ factor in Equation (|1]) set equal 
to 1/2, leading to 

{P{a^ (3)) = ^sm' 26. (3) 

This, to reiterate, can be due to either distance or energy spread or both. 

We now briefly review the various possible solutions to the atmospheric and solar 
neutrino problems before considering the issue of sub-dominant neutrino oscillations. 
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4.3. Atmospheric Neutrinos 

SuperKamiokande detects a 50% deficit of /x-like atmospheric neutrinos coming up 
through the earth ( [Fukuda et al. 1998a,b,(j ). They see no deficit of either upward- or 
downward-going e-like neutrinos. The lower energy, downward-going /x-like events are 
deficient, whereas their high-energy counterparts are not. These data can be explained 
by close-to-maximal z/^ —* oscillations with x ^ e and x = t oi x = s (sterile), or a 
combination thereof. These two alternatives both require parameters in approximately (see, 
e.g., Il'bot, Volkas fc Yasuda 1998|; [I'brnengo, Gonzalez-Garcia, fc Valle 20UUD the same range 



([Scholberg 1999|): 



5m\Ti,^ = IQ-' eV 10'^ and sin^ 29^,^ = 1. (4) 

SuperKamiokande currently favors oscillations to over oscillations to a sterile neutrino 
at the 2a level ( [Fukuda et al. 20"00|) , though some have cast doubt over this result (|Foot| 



2000b| ; pobayashi fc Lim 20001 ). Furthermore, fits to oscillation to a combination of u^- and 



Us fit the data equally well as oscillation to Ur alone ( [Fogli, Lisi, fc Marrone 2001|) . In 
this work, however, for definiteness we will take it that the resolution of the atmospheric 
neutrino anomaly lies in maximal — > u^- oscillations over the relevant length scale, Latm- 

Substituting 5m\^j^j into Equation (|]), we see that the — > oscillation length 
is orders of magnitude less than the size of a typical SNR shell for the entire neutrino 
spectrum (i.e., up to an energy of ~ 1 PeV). This means that the oscillations will be 
distance averaged, and hence at earth we expect a 50/50 mixture of i/^ and Ur from a SNR. 
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4.4. Solar Neutrinos 

The solar neutrino problem can be solved by z/g Vy oscillations, where y = fi,T, s, or 
a combination thereof, are all allowed, with one important proviso: if the Los Alamos LSND 
experiment is correct, then — > oscillations, with parameters that cannot solve the solar 
neutrino problem, have already been detected ( [White 1999|) . So, if the still-controversial 



LSND result is correct, then y = fi is ruled out. The MiniBOONE and BOONE experiments 
at Fermilab should eventually settle this issue (jBazarko 19991). In this work we disregard 



the LSND result, though in the generic six neutrino mass eigenstate model discussed below 
it could be easily accommodated. 

The precise oscillation parameter space required to account for the solar data depends 
on which of the solar neutrino experiments are held to be correct. The two parameter 
ranges defined below, however, are broadly consistent with all solar data; 

1. z/g Uy with a small mixing angle (SMA) 6ey is possible through the MSW effect. 
If this pertains, then the oscillation amplitude will be far too small to affect SNR 
neutrinos, but the most recent data from SuperKamiokande disfavor this solution 
( [Fukuda et al. 20011) . We therefore disregard this possibility. 



2. Ue with a very large mixing angle (LMA) sin^ 29 ey ~ 1 is an interesting possibility 
for the range0 

Sml = 10-3 ^ 10-^° eV\ (5) 



■^Note, however, that some of this range is excluded because of the non-observation of a 
'Day-Night Effect' in the solar neutrino signal, i.e., the absence of the diurnal variation due 
to the MSW-induced 'regeneration' in the earth expected if 6171^ were in the ~ 10^^ eV^ 
range and maximal mixing applies ( |Guth, Randall, fc Serna 1999| ; [Crocker, I'bot, fc Volkas 
T9991; Ibbot 2000a|. 
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The lower end of the solar Sm"^ parameter space, Sm^ <^ 10^^ eV"^ defines 'just-so' 
oscillations where the oscillation length for solar neutrinos is of order 1 A.U. For larger 
StjIq values completely averaged oscillations, with a fiux suppression factor of 1/2 sin^ 
result. Maximal mixing explains almost all of the data with averaged oscillations (excepting 
the Homestake result; [Cleveland et al. 1998|) . Values of 6mQ > 10^^ eV"^ are ruled out by 



the non-observation of Uf, disappearance from reactors (CHOOZ; [Bemporad 1999] and Palo 
Verde experiments; [Boehm 1999| ). 



If the solar problem is resolved with 6171^ <^ lO^^eV"^ range, one finds from Equation 
that Lq becomes larger than an SNR with typical linear dimension ~ 10 pc. within the 
energy range spanned by SNR neutrinos. This means that above some critical energy the z/g 
beam from such a source would be fiavor-coherent (of course, conversely for larger ^rng's, 
SNR would be distance-averaged over their entire energy range). 

In principle, such coherence would evidence itself by an energy dependent spectral 
distortion; the Uf, fiux at a particular energy {E E + AE) would depend on the part of 
the neutrino oscillation wave (for that particular energy) encountered by the earth at its 
distance from the particular SNR. Pragmatically, given the small statistics that will accrue 
from the SNR neutrino sources and the limited energy resolution expected to be achieved by 
any of the proposed neutrino telescopes, one expects no observational consequence of this 
fiavor coherence. This is because the energy dependence of the fiux suppression washes out 
with the inevitably large size of the energy bins particular neutrino events are accumulated 
into. The beam, therefore, is indistinguishable from one in the distance- averaged oscillation 
regime. 

Again with reference to Equation (^, note also that at the extreme lower end of 
the allowable parameter space, StUq = 10"^° eV"^, for a close-by SNR at 1 kpc and at a 
neutrino energy of 1 PeV, the z/g Vy oscillation length, Lq, reaches the order of the 
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SNR-earth distance, Lsnr- The hypothetical SNR's flux at the earth would, then, rise 
from being suppressed somewhat below a PeV to unsupressed somewhat above a PeV if 
the z/g's attained this energy. Unfortunately, we do not reasonably expect any observational 
consequence from this potentially interesting phenomenon; in the flrst place, we would be 
relying on an unlikely coincidence of extreme energy and absolutely minimal allowable 6m?, 
and secondly, neutrino telescopes almost certainly will not be able to detect u^s of SNR 
energy. 

This z/e-bliiidness basically stems from the fact that electrons resulting from charged 
current Ue events of PeV energy and below are arrested extremely quickly in the detector 
medium. Such z/g's are also well beneath the Glashow resonance (|Glashow 1960|; [Berezins"^ 



|fc Gazizov 1977| ; [Gandhi et al. 1996,1998| ) and also the energy thresholds of 'alternative' 



astrophysical neutrino detection techniques that may be suitable for ultra high-energy 
astrophysical detection like air shower arrays ( [Gapelle et al. 1998| ) and radio detection 



of neutrino interactions in Antarctic ice ( paisser, Halzen fc Stanev 1995| ; [Alvarez-Muniz. 



Vazquez fc Zas 1999| ; [Alvarez-Muniz fc Zas 1999]) . The short electron path length leaves 



the experimentalist with essentially no directional information to distinguish signal and 
background. Furthermore, even if directional information could be somehow obtained, the 
electromagnetic showers generated by braking electrons and the accompanying hadronic 
showers from the charge current interactions with nuclei in the detector medium are 
essentially indistinguishable from the purely hadronic showers resulting from neutral current 
interactions of all u flavor at these energies. They are also difficult to distinguish from 
lower energy {E^^ < 100 TeV) charged current events ([The ANTARES Gollaboration| 



1999). In such events the exiting r does not have a large enough 7 factor for it to decay 
sufficiently far from the interaction vertex for the 'double bang' ( [Learned fc Pakvasa 1995[) 
- the hadronic showers due to the original charge current interaction and the later r decay 



to be resolved (see [Crocker, Melia, fc Volkas 2000| and references therein for more detail 
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on the chances for SNR-energy z/g detection). 

Given that double-bang detection of z/^ also only becomes practicable in the uppermost 
decade of the SNR neutrino energy spectrum ([I'he AJNTARES Collaboration 1999| ), in this 
work we will only consider possibilities for SNR detection, though note that SNR u^'s 
remain of some importance because of what they might oscillate into. 



4.5. Atmospheric and Solar Data Combined 

Note firstly that any neutrino mixing scheme that seeks to accommodate both the 
atmospheric and solar neutrino anomalies, with their different attendant oscillation length 
scales, and, therefore, different (5m^'s, must possess at least three distinct neutrino mass 
eigenstates. Furthermore, with just the minimal three neutrino mass eigenstates, if we 
demand maximal ly^^ Vt oscillations over the atmospheric scale and maximal mixing also 
in the solar case, we find that z/g is forced to be maximally mixed with both and Vr over 



solar length scales. In other words, we have bi-maximal mixing ([Vissani 199?] ; [Barger et al 



19981 IBaltz, Goldhaber fc Goldhaber 1998|; IJezabek fc Sumino 19981; lAltarelli fc Ferugliol 



TUggt IMohapatra fc NussinovT^ ). 



That close-to-maximal mixing is demanded by the atmospheric neutrino anomaly 
( [b'ukuda et al. 1998a, b,c ; ApoUonio et al. 19981 ) and also favored by the most recent solar 



neutrino data ( [I'ukuda et al. 2001| ) raises hope for extracting interesting particle physics 
from astrophysical neutrino phenomenology. This is because the poor statistics of the 
proposed neutrino telescopes mean that only modes with large mixing angles, 9, can in 
practice be probed. 
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4.6. Finding Oscillations 

Generically, neutrino oscillations may be evidenced in three ways, viz: 

1. By a difference between the neutrino flavor ratios at point of generation to those 
determined at point of detection, i.e., in the case of SNRs, by a difference between 
i^e : : = 1 ■ 2 : at source and : : Ut — 1 x P(fe I'e) + 2 x ^(1^^ — > 

Ue) : IX P{Ve ^ l^u,) + 2 X P{iy^ ^u) ■ 1 X P{Ve ^^t) + 2 X P{Vf, Vr) at 

detector. In practice, this could probably only be probed by measurement of $°^'^/ ^"^^ 
(given the difficulty with astrophysical detection outlined already), which in the 
absence of oscillations should be zero. 

2. Via an observed discrepancy between detected and expected flux, given one has an 
accurate flx on the absolute flux of a particular neutrino species (which, as already 
discussed, could not be Vf. in the SNR case) expected on the basis of an SNR's 7-ray 
signal: (^"^^ j ^ 1^ where 'l>°^* denotes the observed flux of SNR z/^'s or z/^'s 
and ^^'^^"^ denotes the flux of either flavor expected in the absence of oscillations. 

Essentially, the observation of u^s from an SNR constitutes evidence for oscillations 
under either oscillation diagnostic (1) or (2). 

3. Through the observation of a spectral anomaly, i.e., a distortion of a particular 
neutrino flavor's energy distribution away from its expected shape. In the case of an 
SNR's V signal, this would mean a deviation away from pure power-law scaling of 
neutrino flux with energy (given that it could be determined a priori that the region 
of the distribution under investigation should be governed by a single spectral index). 
The great advantages of this third method over the former two are that: 

i) it only requires observation of a single neutrino species: i/^'s for SNRs in practice. 
One is not obliged to positively identify the v flavor to which the i/^u's might oscillate 



to uncover oscillation evidence nor is one required to know precisely what the expected 
flux of z/^'s is; 

ii) it can give a range for the Sm"^ governing the oscillations, not just a lower bound. 



4.7. Sub-dominant Oscillations 

Now, we have already seen that for the entire allowable Sm\rpj^j and SttIq regimes we 
pragmatically expect totally averaged oscillations in SNR signals, i.e., precisely half the 
z/e's and z/^'s generated in a SNR should oscillate to something else on their journey to 
the earth given bi-maximal mixing (though the actual number of u^s detected would not 
vary from the naive expectation because of oscillations from the maximally mixed Vpr^r 
sub-system to z/g- see, e.g., |Bento, Keranen, fc Maalampi 2000| ; [Athar, Jezabek, fc Yasuda 



2000 ). In other words, the three oscillation diagnostics outlined above, when applied to 
SNR neutrino signals, would only ever act as confirmatory to existing terrestrial neutrino 
oscillation experiments, albeit over vastly different energy and distance regimes. In more 
detail, with the three mass-eigenstate bi-maximal mixing outlined above, we expect (i)z/e : 
z/^ : z/^ = 1 : 1 : 1 at detector with the (possibly) practical, phenomenological implication 
that ^f^^ ~ 1, (ii) (^"J"^ / (^^^^"^ = 1/2 independent of energy, and (iii) no spectral 
distortion of the signal (i.e., a plot of neutrino flux versus energy on a log- log scale 
simply produces a line of constant slope). 

Note, however, from Equation (^, that, if an extra Sm^ significantly smaller than 
10^^° eV^ were operating in astrophysical neutrino oscillations - thereby introducing an 
extra, longer oscillation length scale - things become interesting. In particular, an energy 
dependent spectral distortion might be observed in the SNR's z/^ signal - see later. We label 
this new 5m^, 6m'gjj^, where 'SUB' denotes sub-dominant oscillations. 
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Of course, any scenario invoking this new oscillation scale demands one or more new 
neutrino mass eigenstates. This results in the presence in the theory of additional weak 
eigenstates which are obliged to be sterile neutrinos. Of course, these new weak eigenstates 
present us with alternative possibilities for resolving one or other of the existing neutrino 
anomalies, e.g., we might solve the solar anomaly with fg— > Vg oscillations. 

Consider momentarily the range of Smgu^ potentially probed by SNR neutrino signal. 
To see the spectral distortion mentioned above we require that it occurs, for a particular 
SNR at some distance Lsnr, at an energy below the maximum attained by the SNR 
neutrinos (in fact, somewhat below given the expected tail-off of high energy events with 
a power-law spectrum) and above the energy at which the SNR's signal becomes invisible 
because of atmospheric neutrino background. As mentioned previously, our calculations 
show that this energy is around a TeV for all nine SNRs so far discussed. The Smgu^ 
ranges probed, with such a threshold energy, are approximately 

10-^ ^ 10-^W for 1 kpc and 10"^° ^ lO'^W for 10 kpc. (6) 

One may note immediately that the Suigu^ ranges discussed, given they are so tiny, 
are not ruled out by any existing neutrino oscillation experiment (even with the largest 
5m|[/g, at an energy of 1 MeV, Lsub is still over a million kilometers). Put another way, if 
such Smgjjs scales operate in nature, we can only probe them with astrophysical neutrinos. 
These values for Srrig^^ are tiny numbers. We now go on to discuss how we might motivate 
them. 



- 19 - 



4.7.1. 6-Neutrmo Maximal Mixing Schemes 



As in the case of solar and atmospheric scale oscillations, for there to be an observational 
consequence of sub-dominant oscillations for astrophysical neutrinos, we require that the 
mixing behavior described immediately above be close-to-maximal. We therefore seek 
neutrino mixing scenarios that not only naturally incorporate tiny mass splittings, but also 
result in maximal or close-to-maximal mixing. There are, indeed, two examples of such that 
we know of: the Exact Parity Model ( ll-bot. Lew fc Volkas [boot. Lew fc Volkas 199^ 



[boot 199^ ; [boot fc Volkas 19"^5D and the generic, active ^ sterile, pseudo-Dirac scenario 
( [Wolfenstein 1981 ; Bilenkii fc Pontecorvo 1983 ; Bilenkii fc Petcov 1987 [Kobayashi, Lim ^ 



[Nojiri 1991[ ; [Giunti, Kim, fc Lee 199'^ ; [Bowes fc Volkas 1998[ ; [Geiser 199^ ; [Kobayashi jS[ 
Lim 2000[) . In both these scenarios pairwise maximal and close-to-maximal (respectively) 
mixing between every active neutrino flavor and its sterile partner can naturally explain 
the solar and atmospheric anomalies with and z/^— > . Alternatively, both can 

also accommodate close-to-maximal intergenerational mixing so that, for instance, the 
atmospheric anomaly be resolved by t'^— > Vt oscillations ( [Voon fc Foot 2000 ; Kobayashi fc 
Lim 2000 ), as we assume in this work, whereas the solar neutrino anomaly continues to be 



explained by z/g— > vj , i.e., active to sterile, oscillations. In this latter situation, demanding 
that the scale of the mass-splitting between z/^ and z/^' does not interfere with atmospheric 
neutrino experiment results only constrains ^fn^suB be somewhat less than 10~^ eV"^. It 
is natural, however, to assume that the mass splitting between z/^ and u^' , 5rn?g^Q, be of the 
order of that between z/g and uj , StjIq. 



4.7.2. y^-A/^ewMno Maximal Mixing Schemes 

We are not aware of any particularly strong theoretical motivations for maximal mixing 
in a four neutrino system with very small mass-splitting, but we consider this generic case 
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for completeness. Introducing a fourth, light neutrino mass eigenstate with a mass very 
close to one of the existing states will result in sub-dominant oscillations of both and Ve 
to a new Ug over long length scales. Again if this mass difference is in the Sm^u^ range 
described above and the mixing is maximal, there will be phenomenological consequences 
for astrophysical neutrino signals from SNRs. 

With three mass splittings, 5m\rpj^^ SuIq, and SntguQ and four mass eigenstates, 
there are six possible arrangements of the latter. These can be broken down into two 
'double-doublets' (arrangements of two pairs of mass eigenstates defining SuIq and 5mguQ 
split by Sm\j^j^) and four 'mixed' arrangements. 

5. Observational Consequences of Oscillation Scenarios 

5.1. In Theory 

The phenomenological consequences of the maximal, sub-dominant four and six mass 
eigenstate scenarios mentioned can be gauged by noting the energy dependence that they 
introduce to F^, , the fraction of the total initial neutrino flux from an SNR that arrives at 
earth with flavor i/^: 

F,^ = l/3xP{u,^iy^) + 2/3xP{u^^iy^) = p - asin^ A{E)sub, (7) 

where we define Ascale = ^'^''scale^Lsnr SCALE e {ATM, Q,SUB}. p is 1/3 for 
all the sub-dominant oscillation scenarios under investigation whereas o" is 1/3 in the case 
of the six mass eigenstate scenario, 1/6 for the two, four mass eigenstate double-doublet 
scenarios, and 1/12 in the case of the four, four mass eigenstate mixed scenarios. 

The above expression for F^^ is, in principle, dependent on all three oscillation scales 
- Asub,Aq, and Aatm- In practice, no dependence of F^,^ on Aatm is evident because. 
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over the entire energy range of any SNR, atmospheric oscillations will be averaged as 
already discussed. More interestingly, any dependence on Aq actually cancels out between 
1/3 X P(z/e — > z/^) and 2/3 x P(z^^ v^) essentially because z/g and are maximally mixed 
over the solar and atmospheric scales respectively. 

The energy dependence in p — asin^ A{E)sub will show up in the latter two oscillation 
diagnostics provided sufficient statistics can be accrued and ^to^sub ^^^^ within the range 
defined by (H): over the energy range of an SNR's detected v spectrum, /^^^^"^ will 
go from some constant fraction (< 1/2) well below -Ecrit, to exhibiting oscillatory behavior 
around E = Ecru, to a constant value of 1/2 well above Ecru- Also, in these three 
regimes, (over increasing energy) a plot of the SNR's differential u fiux versus energy on a 
log-log scale produces a line which at first has some constant slope, a, goes through some 
oscillatory regime around E = E^^t and then resumes along the initial constant slope, a. 
By Ecrit we denote the energy at which LsuB{Ecrit) = Lsnr, where Lsub{E) = t^-^c 

"'"■SUB 

and we have employed Equation (|l]) setting 6* = |, x = Lsnr & Lose = Lsub- 

Note that because a is constant within each categorization (six mass eigenstate and 
four eigenstate, either 'double-doublet' or 'mixed'), different examples within each of these 
categorizations produce the same phenomenological consequences in terms of the latter two 
oscillation diagnostics. 

The first oscillation diagnostic does not identify sub-dominant oscillations to a z/g, 
however, because with and Vr already mixed over the atmospheric scale, we simply 
expect $°^^/$°'"' = 1/2, independent of energy. 
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5.2. In Practice 

Of course, neutrino telescopes do/will not measure neutrino fluxes directly. Rather 
they seek to pick out - from the Cerenkov radiation they emit - the limited number of 
charged leptons that cross their instrumented volumes of ice or water that are due to 
charged current interactions of astrophysical neutrino near or in these detector volumes. 
This signal must be extracted from the massive background of down-going, atmospheric 
muons that also cross the detector volume. 

In order to show that the techniques we have described above can be useful diagnostics 
for sub-dominant oscillations in SNR v signals we must determine whether the actual 
event rates in a km^ detector due to promising SNR v sources are large enough that plots 
of neutrino flux versus energy for these sources have small enough statistical errors that 
deviations from pure power- law scaling (i.e., deviations from linearity on a log-log scale) due 
to sub-dominant oscillations might be positively identified against unavoidable statistical 
fluctuations. The derivation of fluxes from event rates also requires that a detector's 
energy-dependent response function be sufficiently well characterized and we assume this is 
the case. 

There are plans afoot for the construction of km^-scale neutrino telescopes in both the 
Antarctic ice (AMANDA and its planned extension IceCube) and the deep Mediterranean 
sea (the ANTARES and NESTOR projects). For recent reviews of the status of neutrino 
telescopy and descriptions of the various neutrino telescope proposals the reader is directed 
to Spiering (2001) for the IceCube project. The ANTARES CoUaboration (1999), and Botai 
et al. (2000) for NESTOR. 
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6. Simple Modeling of Detector Operation 

The probability that a high energy muon-type neutrino is detected in a fcm^-scale 
neutrino telescope depends on two factors, viz; approximately inversely on the interaction 
length of the neutrino (Xint) at that energy (which, in turn, depends on the charged 
current cross-section) and approximately directly on the radiation length of the muon (-R^) 
produced in the interaction (|Halzen 1998|) . (We assume here that the linear dimension of 
the detector is small on the scale of i?^.) We can make a rough estimate of the effect of 
these factors by writing down a detection probability multiplier which goes as some power 
of the energy: 

p 

P ^ — ^ ~ AE " 

In our model we employ the values for n and A given by Halzen (1998). 

We have written a FORTRAN code which, from a particular SNR's 10 GeV 7-ray 
flux, extracts event rates in a neutrino telescope of 1 km^ volume, assuming that the 7-ray 
emission is hadronic in origin. The code takes into account earth shadowing effects relevant 
to a particular SNR given its declination and the detector's latitude. These we estimate 
on the basis of the work of Naumov and Perrone (1999). The event rates are given per 
decade energy bin from 10^ — 10^ GeV to 10^ — 10^ GeV. The code determines these event 
rates given various values for 6111^^^^ and for the minimum and maximum determined 
distances to each SNR. The code does not assume that the neutrino signal has to be 
up-coming for it to be detected. Rather, for nadir angle bins from 10° to 170 ^ 180° 
it compares the signal (more precisely, the differential neutrino flux due to the SNR) to the 
atmospheric background over an assumed detector resolution. This background has two 
components: atmospheric z/^'s and, for nadir angles greater than 90°, atmospheric /i's. We 
employ the zenith angle-dependent parameterization of the sea-level v^j, and ^ flux given by 
Lipari (1993) and also employ the results of Lipari and Stanev (1991) to account for the 
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attenuation of muons with their propagation through the Earth. In the muon-attenuation 
sub-routine we have made the approximation that each detector is located at a single, 
well-defined depth below the Earth's surface (1.6 km in the case of a South Pole and 2.4 km 
in the case of a Mediterranean detector). The code starts to record events in a particular, 
10° angle bin (of which, of course, there are 18 within each energy bin) only when the 
signal rises above the background. Because the atmospheric and /x backgrounds go with 
spectral index a = 3.7 in the energy ranges of concern, whereas the sources go with index 
a ~ 2, the background quickly drops away from the signal once it has been surpassed. 

We estimate the detector angular resolution by employing the parameterization 
suggested by the ANTARES Collaboration (1999): 

^ = (Wr ^ 

This function estimates the importance of the three factors that limit the determination of 
the primary neutrino's direction of travel. These are the uncertainty in the angle between 
the incoming z/^ and the resulting fi, the deviation of the fi away from its original direction 
of travel due to multiple scattering and, lastly, the detector's intrinsic angular resolution as 
determined by uncertainties in its exact geometry, etc. The ANTARES collaboration (1999) 
has determined from Monte Carlo simulations that below 10 TeV total angular resolution is 
limited by the unavoidable angular distribution of the neutrino interactions whereas above 
100 TeV it is limited by detector effects. 

The AMANDA project (which will hopefully evolve into IceCube) has to contend with 
the short scattering length of the Cerenkov light in ice, 24 m , as compared to sea water 
at greater than 200 m. Despite this, IceCube will achieve an angular resolution less than 
one degree and perhaps as low as 0.4° ( |Halzen 2000D and we, perhaps optimistically, adopt 



the same parameterization of detector resolution - Equation (^ - independent of detector 
medium (i.e., ice or water). 
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The binning of event rate data into decade energy bins is forced upon us by two factors: 
the expected hmits to the energy resolution of the proposed detectors and low statistics. 
Again, accurate determination of the energy possessed by a muon neutrino (which produces 
a muon observed by a detector) is limited by three factors: uncertainty in the fraction 
of the neutrino's total energy imparted to the muon, ignorance of the energy loss by the 
muon outside the instrumented volume and, finally, the intrinsic energy resolution of the 
detector apparatus itself ( [I'he AlN'l'AKES Collaboration 19991 ). Given these three factors, 
the ANTARES collaboration (1999) has judged on the basis of MC simulations of their 
detector array that they can gauge a muon neutrino's energy to within a factor of three for 
> 1 TeV, so binning simulated data into decades of energy is reasonable. 



7. Results 

Because of their locations in southern skies, the three SNR-like objects that we have 
determined should produce detectable neutrino fluxes at the earth - the Sgr A East, SN 
1006 and SNR RXJ1713. 7-3946 - must exceed the atmospheric muon background before 
they can be visible to IceCube. This means that, whereas all three emerge above background 
at or below ~ TeV for a detector at Mediterranean latitudes (where they are shielded from 
muons for at least part of the day), at the South Pole their signals only emerge above ~ 10 
TeV. The ten-year, no-oscillation event rates in a km^ detector, per decade energy bin we 
determine for these three objects are presented in Table (1). 

These figures should be halved to incorporate the effect of averaged i/^— Ut oscillations. We 
have taken it that the maximum energy attained by a SNR z/ is 1 PeV (~ 1/12 x 10^^ eV). 

Note that, above background and below the maximum neutrino energy, there are only 
two decade energy bins operating in a South Polar detector for all three neutrino sources. 
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detector 
location 


object 


10^ 10^ 


10^ 10^ 


10^ 10^ 


IqS 10^ 
(in GeV) 


Mediterranean 














Sgr A East 


108 


491 


277 


118 




SN1006 





18 


16 


8 




SNR RX... 





23 


20 


10 


South Pole 














Sgr A East 








156 


152 




SN1006 








1 


11 




SNR RX... 








3 


14 



Table 1: Event rates in decade energy bins due to the three sources under consideration 
assuming no oscillations. 

IceCube might provide evidence for sub-dominant oscillations only via diagnostic (2) - 
observation of energy dependence in the value of $°^^/$^^'^'"'. We repeat, however, that 
positive identification of such variation requires that be well determined which, in 

turn, necessitates accurate determinations of the normalizing differential photon flux at 10 
GeV and the photon spectral index. 

On the other hand, Sgr A East shows up in four energy bins for a Mediterranean 
detector and the other two sources in three. Thus the possibility that the third oscillation 
diagnostic might be brought into play is held out. Whether this diagnostic can be made to 
work in practice depends on statistical error, granted that potential confounds like energy 
dependence in the detector response are well-enough pinned down. 

Figure (2) illustrates the expected, 10-year event rates in a km^ detector at the 
Mediterranean for the most promising oscillation scenario - six mass eigenstate mixing 
- and the best astrophysical neutrino source, Sgr A East. The figure details event rates 
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in decade energy bins with no oscillations and also (atmospheric and) sub-dominant 
oscillations governed by 5m|^^'s from lO"-*^^ to 10~^ eV"^. The error bars indicate the 
statistical error [^/n/n) in each 'datum'. 

Prom Figure (2) one may conclude that ^ml^/^'s larger than ~ 10~^^ eV'^ are not discernible 
one from another. At Srngjj^ values of 10~^ oscillation are completely averaged over Sgr A 
East's observable neutrino spectrum. 

Figure (3) illustrates a way to test the hypothesis that a neutrino spectrum obeys 
a power law. The £^^-weighted differential flux at each particular energy E is related to 
the integral flux from E to 10£^, N[E^ lO^'), one would infer from the event rate (in that 
particular decade energy bin) by the relation 



The figure shows the differential flux at the minimum energy of each of the four decade 
energy bins one might extrapolate from actual observation over ten years of Sgr A East 
with a km^ neutrino telescope at Mediterranean latitude, assuming a power-law spectrum. 
The error bars are calculated from the relative statistical error in the event rate for the 
same energy bin. //the plotted points form a straight line within statistical error, then the 
power law hypothesis is consistent with the 'data'. From Figure (3), however, one may see 
by eye immediately that for a range of Srugu^ values, a constant power- law (i.e., straight 
line) fit to the differential flux is not possible within the error bars. 

Employing a least-squares fit to the differential flux data for the various values of 
SmsjjQ with a two- free-parameter function, f{E) — ax E~^ we determine that in every 
case the fitting algorithm settles to a value for b less than 2. Such a flat spectrum is 
never observed in nature nor is it predicted by shock acceleration models. In fact, a power 
law going with index -2 describes the flattest reasonable spectrum. We therefore fit to a 
one-free-parameter function, f{E) — ax The values determined by such a fit are 



dN{E) 
dE 



1-a N{10E,E) 
IQi-" - 1 E 



(9) 
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presented in Table (2). 



SmQTTR feVl 


d.o.f. 


10-11 


~ 


10-12 


3.2 


10-13 


16.1 


10-14 


19.3 


10-15 


12.0 


no osc 


~ 



Table 2: values for the fitting of a power law of index 2.0 to modelled differential flux 
determinations for different values of drriguQ. 

Note that this oscillation diagnostic works fairly independently of the overall 
normalization of the flux, i.e., one simply needs enough events from the neutrino source 
without having to know exactly its flux in the absence of oscillations. As the table and figure 
above illustrate, Sgr A East does, indeed, produce sufficient events over ten years for this 
method. This cannot, unfortunately, be said for the other two SNRs under consideration. 
Figure (4) shows the differential neutrino flux from SN 1006 for the three active angle bins 
in a Mediterranean km^ detector (generated in the same fashion as for Figure(3) above). 

Again by eye, one may note from Figure (4) that the unavoidable statistical error (calculated 
for ten years' events) washes out the oscillation signature. In fact, the no oscillation case 
may only just be distinguished from all the separate oscillation examples. The neutrino 
signal from SN 1006 (and also that from SNR RXJ1713. 7-3946 and similar, yet undiscovered 
objects), therefore, requires either or both longer observation periods than a decade or 
larger than 1 km^ detectors to be useful in providing evidence for spectral distortion due to 
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sub-dominant oscillations. It may be that the signals from a number of such similar SNRs 
are able to be statistically combined to provide a data set large enough to be probed for 
sub-dominant oscillations. Further, observation (or, indeed, non-observation) of neutrino 
signals from these objects will settle once and for all the question of whether it is hadronic 
or leptonic acceleration that is ultimately responsible for their high energy 7-ray signals 
and will, therefore, be of crucial import to cosmic ray research. 

8. Conclusion 

We have described a general technique via which the effects of sub-dominant neutrino 
oscillations might be uncovered in the spectra of galactic SNRs. This technique does not 
require either observation of i/gS or i/t-'s or unrealistically constrained measurements of these 
objects' high energy 7-ray signals. We have determined through careful modeling that, 
when applied to the Sgr A East neutrino signal, the technique will allow for the discovery 
or exclusion of sub-dominant neutrino oscillations governed by a 5m^ parameter with a 
value in the range 10~^^ to 10~^^ eV^. Such values cannot be probed by any conceivable 
terrestrial and solar system neutrino experiments. 
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Fig. 1. — 7-ray differential fluxes for 7-Cygni and SN 1006. The EGRET data for 7-Cygni 
are from Esposito et al. (1996), the Whipple upper limits from Buckley et al (1998) and 
those due to HEGRA from Prosch et al. (1995). The SN 1006 data are from Tanimori et 
al. (1998) for the CANGAROO collaboration and the EGRET upper limit is quoted by 
Aharonian and Atoyan (1999). The curve passing through the SN 1006 GANGAROO data 
points is the expected integral flux from this object assuming a spectral index of 2.0. This 
curve passes well below the EGRET upper limit. 
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Fig. 2. — Event rates in decade energy bins due to the Sgr A East u source with a number 
of different oscillation scenarios. The horiontal error bars indicate the range of energy which 
defines each of the four active bins. 
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Fig. 3. — Energy^-weighted differential flux from the Sgr A East u source with a number 
of different oscillation scenarios. See the text for a full description of the method used to 
generate this figure. 
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oscillation scenarios. The figure is generated in the same fashion as Figure (3); see the text 
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